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generalized form factors in full QCD with Nf = 2 0(a)-improved Wilson fermions based on simulations closer to 
the chiral and continuum limit. 



1. INTRODUCTION 

Understanding the internal structure of hadrons 
in terms of quarks and gluons (partons), and in 
particular how these provide the binding and spin 
of the nucleon, is one of the outstanding problems 
in particle physics. 

Matrix elements of the light cone operator 
0(a) = /^^g(-^n)^ (1) 

T) —ig f X { 2 , dan-A(an) / A , 

measured in deep-inelastic scattering experiments 
provide a wealth of information regarding the 
quark and gluon content of the nucleon. 

Expanding 0(x) in terms of local operators 
via the operator product expansion generates the 
tower of twist-2 operators 

p^-Mn} =g j 7 {^ D ... D q: ( 2 ) 
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where D= h(D — D) and {• • ■ } indicates sym- 
metrization of indices and removal of traces. The 
(non-)forward matrix elements of Eq. J2J specify 
the (n— l) th moments of the (generalized) parton 
distributions. 

Parton distributions measure the probability 
I'i/'la;)! 2 of finding a parton with fractional longi- 
tudinal momentum x in the fast moving nucleon 
at a given transverse resolution 1/Q. General- 
ized parton distributions (GPDs) |1I2| describe 
the coherence of two different hadron wave func- 
tions ijp{x + £/2) ip(x — C/2), one where the par- 
ton carries fractional momentum x + ^/2 and one 
where this fraction is x — £/2, from which further 
information on the transverse distribution of par- 
tons can be drawn |3I4| . In the limit where the 
momentum transfer A to the nucleon is purely 
transverse, i.e. A = (0, Aj_) and £ = 0, GPDs 
regain a probabilistic interpretation 4 . In the 
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forward limit (A = 0), these distributions reduce 
to the Feynman parton distributions. 

Moments of GPDs are amenable to lattice cal- 
culations |5I6I7| . Thus, they offer a promising- 
way to link phenomenological observations to first 
principle theoretical considerations. In this talk 
we report on recent unquenched results obtained 
by the QCDSF collaboration. 

2. SIMULATION DETAILS 

We simulate with N{ = 2 dynamical config- 
urations generated with Wilson glue and non- 
perturbatively 0(a) improved Wilson fermions. 
For five different values (3 = 5.20, 5.25, 5.26, 
5.29, 5.40 and up to three different kappa values 
per beta we have in collaboration with UKQCD 
generated 0(2000 — 8000) trajectories. Lattice 
spacings and spatial volumes vary between 0.075- 
0.f23 fm and (f .5-2.2 fm) 3 respectively. 

Correlation functions are calculated on config- 
urations taken at a distance of 5-10 trajecto- 
ries using 8-4 different locations of the fermion 
source. We use binning to obtain an effective dis- 
tance of 20 trajectories. The size of the bins has 
little effect on the error, which indicates auto- 
correlations are small. This work improves on 
previous calculations by adding one more sink 
momentum, p 2 , and polarization, T\. We use 
Po = (0,0,0), fx = (;p,0,0), p 2 = (0,p,0) 
(p = 2tt/L s ) and r unpol = ~(1 + 74), Ti = 

5(1 +74)47571, r 2 = 1(1 +74)17572- 

3. QUARK DISTRIBUTIONS 

The moments of the F 2 structure function 
are determined by calculating the matrix 
elements 



(N(p)\ 



\N(p)y 



2v^ s (g s (M)) [^•••p^-Tr], (3) 

and corresponding Wilson coefficients in a renor- 
malization scheme S (eg. MS) and at a scale Ai 
(e.g. 2 GeV). See [8] for the operators to use on 
the lattice. To eliminate 0(a) terms in the matrix 
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Figure f . Bare results for the lowest moment of 
F 2 at (/3,K sca ) = (5.4,0+3560), t = 17. The fat 
horizontal lines are obtained from a fit to the data. 
v 2a an d v 2b refer to different lattice versions of 
the same operator. 

elements, these operators need to be improved. 
For the lowest moment, i.e. v 2 , this amounts to re- 
placing 0^ v by (l+am q c ) O^+EiLi ac i 
The improvement coefficient cq is only known per- 
turbatively, while the other coefficients are un- 
known. Since the corresponding operator ma- 
trix elements turn out to be negligible we will 
set ci = C2 = 0. For the higher moments an ad- 
ditional problem arises: the operators may mix 
with operators with total derivatives [§]. How- 
ever, the corresponding operators are again found 
to be small and are therefore neglected here. 

Matrix elements are determined from the ratio 
of three-point to two-point correlation functions 



Tl{t,T-p',p-0) 



C r (t,r;p',p,0) 



(4) 



C 2 (t,p') 

C 2 (T,pl)C 2 (t,p')C 2 (t-T,p) 
C 2 (T,p)C 2 (t,p)C 2 (t-T,p') 



where C 2 is the unpolarized baryon two-point 
function with a source at time and sink at time 
t, while the unpolarized three-point function Cr 
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Figure 2. Finite size analysis for (p 
(0,0,0); at (/3,/c S ea) = (5.29,0+3550) on V 
24 3 x 48, 16 3 x 32, f2 3 x 32 lattices. 



3 



0.4 
0.2 



L 






0.4 








- m 










36 






V2m K 






(u-d) RGI - 


0.2 

! 





1 


2 3 


4 

(r m ps ) 2 


5 


6 


7 






0.4 


















13 
















0.2 



( 


> 




0.02 




0.04 




0.06 



(a/r ) 



Figure 3. Results from a fit to Eq. |||) with 
The upper plot shows the chiral extra- 
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polation (ro, mps) — C2b(a/ro) and the lower 
plot the continuum extrapolation u^p^ro, mps) — 
F v *»(r m FS ). 

has an operator O insertion at time r. For the 
matrix elements v n we have p = p ', i.e. the last 
term of Eq. is equal to 1 . To improve our sig- 
nal for non-zero momentum we average over the 
results for Cr and C2 for both momenta before 
we calculate the ratio 1Z. In Fig. ^ we compare 
the results for the different ways of obtaining the 
lowest moment. Within errors the results are con- 
sistent. 

The bare matrix elements must be renor- 
malized. While non-perturbatively determined 
renormalization factors are known in quenched 
QCD JOj) the situation for dynamical sim- 
ulations is still under investigation. Here 
we use tadpole-improved renormalization-group- 
improved boosted perturbation theory jl lj to 
convert our lattice results to RGI. 

Before we examine the quark mass and lattice 
spacing behaviour of our results, we first check 
for finite size effects. In Fig. [5] we present the 
first results of a finite size analysis for V2b- Here 
we plot V2b as a function of inverse spatial lattice 
extent for three different volumes V = 24 3 x 48, 
16 3 x 32, 12 3 x 32 at (f3, rc sea ) = (5.29,0.13550). 
This preliminary analysis reveals that the finite 
size effects for i>2& are small. 

Finally, the discretization effects and quark 
mass dependence need to be investigated. The 
commonly accepted procedure is to first extrap- 
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Figure 4. Same as Fig. but for A x = 1 GeV. 

olate to the continuum limit then to the chiral 
limit. Since the currently available data does 
not allow us to perform both extrapolations sep- 
arately, we make the following ansatz for a simul- 
taneous chiral and continuum extrapolation: 



C(^o,mps) 



F"™ (r mps) 

2 



ro 



+ d n arom I>s . (5) 



The terms proportional to c„ and d n take dis- 
cretization errors oc a 2 and residual effects oc am q 
into account (we fix d n — 0). The first term corre- 
sponds to a chiral extrapolation. To incorporate 
chiral physics into the extrapolation taking into 
account the effects due to the 'pion cloud' sur- 
rounding the nucleon, the ansatz 



,,RGI 



1 - cx In • 



■ a n x 



(0) 



has been proposed J21> where c « 0.663 and 
A x is a free parameter, usually taken to be of 
0(1 GeV). Setting A x = reduces Eq. © to a 
linear extrapolation in (r mps) 2 . 

In Fig. we plot the results for v^ b Gl (r , mps) — 
C2b(a/ro) 2 as calculated from Eq. (© with A x = 
(x 2 /d.o.f = 2) and plotted as a function of mp s 
for the isovector u — d. We see here that al- 
though our data at heavy quark masses agree 
very well with a linear extrapolation, the pre- 
dicted value in the chiral limit is roughly twice 



4 



3 


(u-d) RGI - 

— h» •. 




: 


(u-d) RGI - 

- 





0.02 0.04 0.06 



Figure 5. Continuum extrapolation of V3 and v^. 

the experimental value. A fit using Eq. (JSJ with 
A x = 1 GeV and the result in the chiral limit 
constrained to the experimental value does not 
describe our data equally well (x 2 /d.o.f = 5). 
However, as can be seen from Fig.0]any disagree- 
ment is not significant since all of the curvature of 
the fit occurs in the light quark mass region where 
there is no data. Even at our lightest quark mass 
(mps ~ 500MeV) we are not yet into the region 
where curvature is expected to occur. 

Turning our attention now to discretization ef- 
fects, in the bottom figure of Fig. 01 we plot 
u 26 ( r 0i m Ps) — F V2h (ro^ps) as a function of 
(a/ ro) 2 . Here we observe a very small depen- 
dence on the lattice spacing, indicating that not 
only is our 0(a) improvement program working, 
but also that 0(a 2 ) effects are small. 

For the higher moments we only perform a chi- 
ral extrapolation linear in (romps) 2 - The results 
for V3 and V4 from the combined chiral and con- 
tinuum extrapolation are shown in Fig. 

4. GENERALIZED PARTON DISTRI- 
BUTIONS 

Non-forward matrix elements of the twist- 2 op- 
erators in Eq. yield the moments of general- 
ized parton distributions 

J dxx^HgfaCA*) = H qn (£,A 2 ) 
J dxx 71-1 E q (x,£,A 2 ) = E qn (£,A 2 ) (7) 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
-A 2 [GeV 2 ] 



Figure 6. Generalized form factors A 20 d , -B20 d 
and C2 ~ d , together with a dipole fit. 

where [2] 

H qn (Z,A 2 ) = ^A 9 „, 2j (A 2 )(-20 2 * (8) 

i=0 

+ Mod(7i + l,2)C 9n (A 2 )(-20" 

E qn {^A 2 ) = ^> g „^(A 2 )(-2£) 2 * 

- Mod(n+l,2)C g „(A 2 )(-2£)" 

and the generalized form factors A^ 2i (A 2 ), 
Bl 2i (A 2 ) and C«(A 2 ) for the lowest three mo- 
ments are extracted from the nucleon matrix el- 
ements (p'lO^ 1 •""•»> \p) |2j ■ For the lowest mo- 
ment, Aiq and Biq are just the Dirac and Pauli 
form factors F\ and F% , respectively. We also ob- 
serve that in the forward limit (A 2 = £ = 0), 
the moments of H q reduce to the moments of the 
unpolarized parton distribution A n Q = (a;™ -1 ). 

In order to extract the non-forward matrix ele- 
ments we compute ratios of three- and two-point 
function as in Eq. J3J| with p^p'. 

In Fig. we show, as an example, the general- 
ized form factors A20, -B20 and C2 for the non- 
singlet, u — d, on a 24 3 x 48 lattice at j3 = 5.40 and 
tsoa = K vai = 0.13500 corresponding to a lattice 
spacing, a ro — 6.088 and mps ~ 970 MeV. 

The generalized form factors A20, -B20 are well 
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Figure 7. Generalized form factors A^Q d , 

and A^Q d , together with a dipole fit. All form 

factors have been normalized to unity. 

described by the dipole ansatz 

Am 



An{A 



fl - A 2 /M 2 ) 



(9) 



while Co d is consistent with zero. 



Burkardt @] has shown that the spin- 
independent and spin-dependent generalized par- 
ton distributions H(x, 0, A 2 ) and H(x,0, A 2 ) 
gain a physical interpretation when Fourier trans- 
formed to impact parameter space 



q(x,bj_) = 



d 2 A± 

W) 2 



-ib± -Aj 



H(x,0,-&1), (10) 



(and similar for the polarized Aq(x,b±)) where 
q(x, b±) is the probability of finding a quark with 
longitudinal momentum fraction x and at trans- 
verse position (or impact parameter) b±. 

Burkardt @] also argued that H(x, 0, — Aj_) be- 
comes A^-independent as x — > 1 since, physi- 
cally we expect the transverse size of the nucleon 
to decrease as x increases, i.e. lirn. x _>i q(x, bj_) oc 
S 2 (b±). As a result, we expect the slopes of the 
moments of H(x, 0, — Aj_) in A^ to decrease as 
we proceed to higher moments. JThis is also true 
for the polarized moments of H(x, 0, — Aj_), so 
from Eq. @ with £ = 0, we expect that the 
slopes of the generalized form factors A n o(A 2 ) 
and A n o(A 2 ) should decrease with increasing n. 
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Fi gure 8. Same as Fig. ^ but for A^q 

In Figs. [3 and |H| we show the A 2 -dependence 
of A n0 (A 2 ) and A n0 (A 2 ), respectively for n = 
1,2,3. The form factors have been normalized to 
unity to make a comparison of the slopes easier 
and as in Fig we fit the form factors with a 
dipole form as in Eq. 0. We note here that the 
form factors for the unpolarized moments are well 
separated and that their slopes do indeed decrease 
with increasing n as predicted. For the polarized 
moments, we observe a similar scenario, however 
here the change in slope between the form factors 
is not as large. This is to be compared with the 
results from Ref. [7] which reveal no change in 
slope between the n — 2 and n — 3 polarized 
moments. 

Although fitting the form factors with a dipole 
is purely phenomcnological, it does provide us 
with a useful means to measure the change in 
slope of the form factors by monitoring the ex- 
tracted dipole masses (Mi, M 2 , M 3 ) as we pro- 
ceed to higher moments. We have calculated 
these generalized form factors on a subset of our 
full complement of (/3, k) combinations and have 
extracted the corresponding dipole masses. We 
plot these dipole masses in Fig. as a function 
of rripg. Included for comparison are previous 
quenched results for the first ^21 an d second |3] 
moments. The new unquenched results indicate 
that quenching effects are small for form factors. 

The important feature to note in Fig. is 
the distinct separation between (and increase in 
magnitude of) the dipole masses as we move to 
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Figure 9. Dipole masses for the first three mo- 
ments of H (x, 0, — AjJ as a function ofrrip S , to- 
gether with a linear extrapolation to the physical 
pion mass. 

higher moments (x — > 1). Although the data 
available for M3 is limited, the behaviour of all 
dipole masses appears to be linear with »ip S . 
Consequently, we perform individual linear ex- 
trapolations of the dipole masses M%, M2, M3 to 
the physical pion mass, although the findings of 
Rcf. 14 suggest that the chiral extrapolation of 
the dipole masses of the electromagnetic form fac- 
tors may be non-linear. 

If the dipole behavior Eq. J3J continues to hold 
for the higher moments as well, and if we as- 
sume that the dipole masses continue to grow in 
a Regge-like fashion, we may write 

j x ^ a"" 1 H g (x, 0, A 2 ) = (i _ A ? 2/M 2 )2 , (11) 

with M? = const. + I /a', where const, » —0.5 
GeV 2 and 1/a' « 1.1 GeV 2 . This is sufficient 
to compute H q (x,0, A 2 ) by means of an inverse 
Mellin transform |15| . 

Having done so, the desired probability distri- 
bution of finding a parton of momentum fraction 
x at the impact parameter b± can then be ob- 
tained by the Fourier transform of Eq. I|10(l . 

5. CONCLUSION 

We presented an update on our results for 
(a;™" 1 ) in full QCD. While data closer to the chi- 



ral and continuum limit became available, the lin- 
ear behaviour observed previously |16| still holds. 
Finally, we have shown preliminary results for the 
first 3 moments of the GPDs. Our data confirms 
the expected flattening going to higher moments. 
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